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The study of nanospintronic devices utilizing coherent transport through molecular scale multiply- 
connected geometries in the presence of moderate magnetic fields is presented. It is shown how two 
types of simple devices, spin filters and spin splitters (or Stern-Gerlach devices) may be constructed 
from molecular nanometric rings utilizing the Aharonov-Bohm effect. The current is calculated 
within a single electron approximation and within a many-body master equation approach where 
charging effects are accounted for in the Coulomb Blockade regime. We provide rules and tools to 
develop and analyze efficient spintronic devices based on nanometric interferometers. 



I. INTRODUCTION 

In recent years a great deal of attention has been de- 
voted to the study of useful electronic devices utilizing 
the Aharo nov-Bohm^ (AB) phase in multiply-connected 
geometries.'^I^IlIIISEIl] In partic ular , there has been in- 
terest in spin-sensitive device^SEo] g^j-g ^^j-^g single- 
electron analogue of semiconductor spintronics.^^ Most 
of the research conducted in this direction has focused on 
mesoscopic systems, where the AB flux quantum matches 
weak magnetic fields, inter-electronic dependencies can 
be relatively negligible, and the RashbaP^ and Dressel- 
hauj ^ spi n-orbit coupling or inhomogeneous magnetic 
fields,13I151 provide a large and controllable dependency 
of the electron Hamiltonian on directionally significant 
spin eigenvalues. 10 n is 19 20 23221231 

At the nanometric scale, it has been shown that spin- 
independent AB molecular interferometers may be pos- 
sible at reasonable magnetic fields when the coupling 
of the device to the leads is small.l^llllslllll p^j. 3^^}^ 
small AB interferometers there are striking differences 
in the properties of magnetic versus electric gauges, and 
the magnetic gate becomes advantageous over electrical 
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For example, the current in a multi-terminal 



molecular device can be tuned by changing the polar- 
ity of the magnetic field utilizing its symmetry breaking 
nature. Another example includes fundamental differ- 
ences between magnetic and electric gates with respect 
to inelastic effects.l^ 

Previous studies of nanometeric molecular AB inter- 
ferometers have ignored the spin degree of freedom. In 
such molecular-scale systems, if the likes of atomic spin- 
orbit coupling and magnetic ions are not prevalent ,1221201 
the lack of a Rashba/Dresselhaus field and the difficul- 
ties in locally manipulating the external magnetic field 
at such scales leaves only the much smaller Zeeman term 
to differentiate between spins.'^^' In addition, the small 
size of such devices further complicates matters by in- 
troducing large charging effects.'22l These should become 
especially noticeable when the coupling between the de- 
vice and the leads is weak (as required for nanometeric 
devices operating at reasonable magnetic fields), since 
electrons spend a longer time on the device. Despite 



these inconveniences, one may very well wonder whether 
it is possible to develop molecular AB interferometers as 
single-electron spin-devices such as a spin filter or split- 
ter. This is the question that interests us here. 

In this paper we develop the basic physical founda- 
tion to develop such devices. In Section |TT] we describe 
the types of devices we wish to model and define their 
efficiency in terms of the spin-dependent current. This 
includes a two-terminal spin filter and a three-terminal 
spin splitter. In Section |III| Zeeman splitting is intro- 
duced into the two models developed by Hod et al^^^^^ 
for independent spinless electrons - a simple continuum 
one-dimensional model and a more detailed tight-binding 
model. Within the simple continuum model we pro- 
vide an exact solution for the spin-dependent conduc- 
tivity for the two- and three-terminal devices. The role 
of the different model parameters is studied and com- 
parison between the simple continuum model and the 



tight-binding results are made. In Section IV the tight- 
binding levels are used in a master equation calcula- 
tion^^ 34 35 36 37 38 39 40 41 42, gxamine charging effects 
in the Coulomb blockade regime. From this basic anal- 
ysis some interesting conclusions are drawn regarding 
the properties and limitations of nanospintronic AB de- 
vices, and a few potential avenues for further research are 
pointed out in Section \V\ 



II. DEVICES 

We will examine two types of elementary nanometric 
devices with and without charging effects: a spin filter 
and a spin splitter, as illustrated in Fig. [l] Such devices 
might be built from molecular rings such as the polycyclic 
aromatic hydrocarbon,^^ or from atomic corrals 
(with th e adv antage of controlled structure and electron 
density) I '^^ l '^^ l Such corrals can be constructed atom by 
atom using scanning tunneling microscopy techniques.'^ 
A more exotic example involves 3D nanostructures such 
as torus knots, where a periodicity of more than 2tt can 
be achieved, producing a stronger AB effect than sim- 
ple rings. In such molecular structures, control over the 
structure's geometry is possible at the atomic level (un- 
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FIG. 1: A sketch of the two prototype devices: a molecular 
spin filter (upper panel) and a molecular spin splitter (lower 
panel). Up/down refers to the two spin polarized states. 



like the case of mesoscopic rings) . 

The filter does not really require something as com- 
plex as an AB ring for its realization. One only requires 
two single-spin levels which can, by manipulation of the 
magnetic field, be moved in and out of the conduction 
energy window, and any atom, quantum dot, or other 
discrete system can meet this requirement. This case is 
treated because it provides a simple example with only 
two terminals and yet may still be useful. The splitter is 
a molecular Stern-Gerlach device and is somewhat more 
interesting: having three terminals, it actually requires 
a device which, like an AB ring, has some sort of inher- 
ent asymmetry in either its construction or in the spinor 
wavef unction. 

It is useful at this point to examine the energy scales 
of the problem. The Fermi energy is the most important 
parameter. For a half-filled conduction band the Fermi 
wavelength is of the order of four bond lengths Xp — 4a. 

7 2 

This gives a Fermi energy oi ep — ^m* ~ lOeV if 
we set the lowest level to zero, regardless of ring size and 
assuming atomic spacing of about a — 2 A. To a good ap- 
proximation, the molecular orbital energy (at zero mag- 
netic field _B = 0) is given by an effective mass model 
— 2m*'R^ ' where m = 0, ±1, • • • , i? is the radius of the 
ring and m* is the effective mass. This gives a value for 
the Fermi quantum number of = — . Electronic lev- 
els on the rings near the Fermi energy are separated by 

^2 

Ae = jn*Ra ' '^hich reduces linearly with R. For rings of 
the order of several nanometers in diameter the separa- 
tion between levels near the Fermi energy is of the order 
of sub electron volts. This is also the order of magnitude 
of the charging energies Ec. Other energy scales involve 



the magnetic field splitting which can be as large as the 
level spacing and the Zeeman splitting 2.89 • 10~^eV per 
Tesla, which is only about a third of the thermal energy 
at IK. 

The effectiveness of these devices can be judged not 
only by output polarization, but also by unitarity: a per- 
fect device always sends the correct electrons into the 
correct lead. More precisely, the two devices that will be 
discussed are judged by the following quantities: 
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Here is the up/down current for the two terminal de- 
vice and is the up/down current for channel 1 or 2 for 
the three terminal device. The only important proper- 
ties these quantities need to have for the purpose of this 
work is that they should be bounded from above by 1, 
and should reach 1 only in the case of perfect operation. 
For the two and three terminal devices both requirements 
are satisfied. We refer to them from now on as "efficien- 
cies" . 

The problem of building a perfect device can thus be 
mathematically restated as the problem of optimizing its 
efficiency to the desired value, usually unity, over the 
space of all controllable parameters. These include the 
magnetic field, the ring's radius and structure, the gate 
and bias voltage (or the chemical potential) and the lead 
placement and coupling. Reasonable ranges for these pa- 
rameters must be assumed: magnetic fields of more than 
a few Tesla may pose a technical limitation, as are large 
rings where disorder effects begin to dominate. The 
leads cannot be too close to one another to avoid direct 
tunneling and should be coupled strongly enough to make 
the environmental coupling unimportant. Even within 
these limits, the problem remains numerically formidable 
in all models. We will therefore also discuss a conceptu- 
ally simpler if less systematic way of designing perfect 
devices. 



III. SINGLE ELECTRON PICTURE 

In this section we describe the physical principles re- 
quired to construct molecular spin filter and spin splitter 
devices in the limit where charging effects can be ignored 
and the complexity of the many-body physics can be re- 
duced to a simplified one-electron picture. 



A. Basic Concepts 

To discuss the basic physical principles involved, let us 
first regard a simple analytical model for the transmission 
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through an AB ring! ^° ' ^^ ' ^^ U^'22lThe physics of this model 
for molecular conductors was recently discussed by Hod 
et al. in the absen ce of Zeeman splitting and spin-orbit 
coupling .'22IIII2IIII1 Consider a one dimensional AB inter- 
ferometer as shown in Fig. [l] It consists of a conducting 
ring of radius R coupled to two or three conducting wires 
placed in a perpendicular uniform magnetic field, B. The 
device is described by the Pauli Hamiltonian: 

ffs = ;^^ + nr)+.WB^-B, (2) 

where 11 = P + ^A(r), A(r) is the vector potential (B — 
V X A), V{y) is the electrostatic potential on the ring, 
and TOe is the electron mass. The last term in Eq. ([2| 
represents the coupling of the magnetic field to the spin 
angular momentum, where <t are the Pauli matrices, g is 
the gyromagnetic ratio (we take g = 2 for the spin), and 
/1b = is the electron Bohr magneton. The lack of a 
Rashba iield in molecular conductors allows us to safely 
ignore spin-orbit coupling effects. In the common case 
where the scalar potential V{r) governing the system is 
periodic or nearly so, the effective mass approximation 
can be used and the Pauli Hamiltonian can be reduced 
to: 

Hs — ^ + ■ B = Hring + (3) 

Zm* 

where m* is the effective mass of the electrons. Since the 
two additive terms in the Hamiltonian must commute 
(one depending only on the orbital part and the other on 
the spin part of the wavefunction) , we can adopt the so- 
lution of Hod et al?^ for the transmission of electrons 
through the two prototype devices (spin filter and split- 
ter). We adopt a transformation e — > e — Eo- which allows 
us to project spin effects onto a calculation where spin 
was previously neglected, noting that the spin dependent 
part of the Hamiltonian under a constant field must have 
the two eigenvalues eo- = ±<7/^b^ for spin up/down elec- 
trons. This transformation is not limited to the specific 
Hamiltonian described above: the only requirement is 
that Hring docs not contain a spin dependency. Further- 
more, even for inhomogeneous magnetic fields, where the 
spin-dependent term in the Hamiltonian does not com- 
mute with Hring ■ ^ similar transformation in the limit of 



adiabatic spin dynamics can be made, where in addition 
to a shift in energy one has to introduce a shift in the 
magnetic 

Using the standard analytical approach of treating a 
ID ring based on a scattering matrix formalism,!^ the 
transmission as a function of en ergy, previously calcu- 
lated for spinless electrons,'^^'2ll need only be modified 
by the Zeeman energies for up (down) electrons, which 
affects only the kinetic phase angles (pk in the expression: 

4,1} =(^k{^\{) ^ ^y^2m*e'^^ = '^^2m*ie± g^i^B). 

The transmission itself is the solution of the linear scat- 
tering problems with the final results for the two-terminal 
device shown in the upper panel of Fig. [l] given by: 

Tnie)^—, ^nil + cos24nr) 

-^Ti ^Ti -I- Q^i cos2 2(t>„i 

where 0„i = tt^ is the ratio between the magnetic flux 

(f> = t:R^B and the quantum flux = and we have 
defined: 

A^^ = l&e\l-cos24l^) 

= 2(c- l)2(c+ 1)2 - 4(c2 + l)(c+ l)2cos2,^^^ 

gu-(c+i)4 

i?U = (c-l)4 + 4c4 + 4- 

4(c2 -I- l)(c- l)2cos24^ -l-Sc^ cos44^ 

(6) 

with the transmission amplitude into the junctions 
and c ~ ^/l — 2e the junction scattering amplitude.'^ 

A similar calculation can be made for the more cum- 
bersome case of a three terminal device shown in the 
lower panel of Fig. [T| We focus on the case where all 
three junctions have identical scattering amplitudes. The 
transmittance for channels 1 and 2 is given by the ra- 
tio T^fie) = N^f{e)[D^^^{e)]-\ The denominator of 
the transmittance probability for both output channels 
is given by 
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FIG. 2: Plots of the conduction (in units of the quantum conduction go = ^) of a spin filter device computed within the 
analytical model at zero temperature and zero bias. Left: the conduction as a function of the magnetic flux for different model 
parameters. The parameters used are: R = 5nm, e = m* = me and Xf ~ IsA (solid lines) and \f = 4A (dashed lines) 
for the lower left panel, R = 5nm, e = j^, m* = me and Xp = 12A for the upper left panel, R — lOnm, e = m* = me 
and Xf = 12 A for the lower right panel, R = 5nm, e — j^^, m* — lOme and Xf ~ 12 A for the upper right panel. Right: The 
spinless conductance as a function of the energy and the magnetic flux. The value of the conduction for the spin up/down are 
given by the values at the ends of the white lines. This figure contains the same information as that shown in the left panels. 
Rainbow color scheme used (red indicates high conductance of ~ Igo)- 



D^^{e) = — (c^ + 1)(19- 12c + 2c2 - 120^ + 19c^) + 320^ cos(40^^) 



2(c- l)'^c{cos[44'^(l - 2a)] +cos[4<?!)[^(l - 2/3)] + cos[4(/)[Hl - 27)]}- 



c- l)^c(c^ + l){cos[4<?!>J^J-(a- 1)] +cos[(40|.^/3- 1)] + 



Ti, 



'^1^(7-1)]}- 



4(c - 1)2(2 - c + 2c^ -c^ + 2c^)[cos{A<pl^a) + cos(40^^/3) + cos(44^7)] + 
2(c - l)4(c2 + l){cos[40^i(a _ p)] + cos[A^l^{a - 7)] + cos[40^i(^ _ ^)]|_ 

^(c + 1)* {-4[1 + c(c - 1)] cos(20^i) + (c - l)2{cos[20^^(l - 2a)] + cos[2,^^^(l - 2/3)] 
+ cos[20l^(l - 27)]}}cos(20„) + l(c + l)^ cos2(20™). 



(7) 



The numerator of the transmittance probability through output channel 1 is given by 



iV|^(e) = -^e^ {-4(1 + c2)2(c - 1)^ cos(44^a) + (c + 1)' cos(40^J-/3) + 2(c - 1)^ 008(44-^7)+ 

(a + 7)1 -(c- 1)2 cos [40^^(a-7)l -2c(c+l)cos \2(t)„^ - 2(1)1^] -2(c+l)cos [20,„ + 20^^^ 



4c cos 
(c^ - 1) cos 
2c(c + 1) cos 



2(t),n + 2{l ~2a)(t)l^ +(c2-l)cos 20„ - 2(1 - 2a)0],-^ +2(c+l)cos 



20™ - 2(1 - 2(3)4^ - {c' - 1) cos 2<j)m + 2(1 - 2j)(^l' + [c' ~ 1) 



20™ + 2(1 -2/3)01; 
2(/)m - 2 



(l-27)'/'I>]} 
(8) 
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and the numerator of the transmittance probability through output channel 2 is given by 
Nf^{e) = -^£2 |_4(i + + {c+lfcos{A(Pl^a)+2{c-lfcos{4<pl^f3)+2{c~lfcos{A<j>l^j)+ 



4c cos 



40l^(/3 + 7) -(c-l)^cos 401^^ 7) -2c(c+l)cos 



Ui 



2c(c-|-l)cos 2(?!)„ + 2(1 - 2a)(/)J,-^ +2(c+l)cos 20„ - 2(1 - 2q;)0].-^ +(c^-1)cos 



2(j)jn + 20 



2(c+ 1) cos 



(c^-l)cos 2(/)„ - 2(1 - 2/3)(/)[> +(c^-l)cos 2<?!)„ + 2(1 - 27)0[^ -(c^-l)cos 



2(/)„ 

20m 



20™ - 20j,"^ 
-2(l-2/?)0Ti 
2(l~27)0Ti]} 



(9) 



In the above equations a, /? and 7 = 27r — a — /3 are the 
angles between all three leads as defined in Fig. [T] 

The current is related to the transmittance through 
the Landauer formula.!^ For the two terminal device the 
current is given by 



e 



de{fie-^lI)~fie~^^o)}Ttl{e) (10) 



where /i/,o £^re the chemical potentials of the in- 
put/output channels and the Fermi function is /(e) = 
, , },t T ■ For the three terminal device the current is 
given by 



■1.2 

hi 



e 

■nh 



de- {/(e - M/) /(e - MOi.O, 



(11) 

where ^iOi,02 is the chemical potential for output channel 
1 or 2, respectively. Conductance can be obtained from 
the current by taking the derivative with respect to the 
bias voltage. 



B. Spin Filters 

We now turn to discuss the application of the above re- 
sults to the construction of a spin filter device. While f or 
the case where the Zeeman effect was neg 

lecte(Pl2ll23 it 

was always desirable to vary the gate voltage so as to shift 
the conduction peaks near zero magnetic field, in Zee- 
man spin devices this heuristic is complicated by the de- 
pendence of the splitting on the magnetic field strength. 
In practice, a compromise between realistically low mag- 
netic fields and usable energy shifts restricts the sought 
after set of parameters. These include the kinetic phase 



angle 



and the transmis- 



sion amplitude into the junction ^/e (cf., Eqs. ^ and 
([6])). Note that the kinetic phase angle depends on the 
magnetic phase angle 0m- This is precisely where the 
aforementioned complication enters. 

With the above convenient expressions one can calcu- 
late the spin-dependent and spin-independent conduction 
for different ring configurations and external parameters. 



Typical results are shown in Fig. [2] for zero bias volt- 
age /i/ = fiQ E^nd for T = OK. As expected, the spin- 
independent conduction is periodic with a period that 
is equal to 0o- Therefore, we plot only the first period, 
namely < 0/0o < 1- We find that within each pe- 
riod the spin-independent conduction has a symmetric 
struc ture around 0/0o = 1/2, characterized by a double 
peak.'^^'^^' This structure is caused by resonance trans- 
mission through the energy levels of the ring. The spin- 
dependent conduction follows closely the behavior of the 
total conduction. The up/down conduction peaks are 
separated by the Zeeman splittings which increase with 
increasing magnetic flux. 

Similar to the case of spinless electrons,'^ changing the 
coupling strength between the leads and the ring, mod- 
eled here by e, changes the width of the peaks as shown 
in the left panels of Fig. |2] The effects of changing the 
kinetic phase angle are more involved than the case stud- 
ied before.t^^In the case of spinless electrons, the kinetic 
phase angle depends only on the product R^/mFe. There- 
fore, changing the energy of the conduction electron by 
the application of a gate voltage, changing the ring di- 
mensions, and changing the effective mass of the conduc- 
tion electron can be mapped onto a universal curve.^^S^ZI 
In the present study there are three independent param- 
eters that affect the kinetic phase angle in different ways. 
For example, the position of the conduction peaks can be 
shifted by adjusting the energy of the conduction electron 
by the application of a gate voltage. This is depicted in 
the lower left panel of Fig. [2] However, the change in the 
energy of the conduction electron also affects the splitting 
between the conduction peaks of the up/down spins as 
previously discussed. Altering the effective mass within 
the ring with the chemical potential held constant, for 
instance by a change of composition or inter-atomic dis- 
tance in the ring atoms, modifies the kinetic phase and 
can thus increase or decrease the splitting as shown in the 
upper right panel of Fig. [2] Finally, as the ring radius is 
increased, the splitting for a set flux ratio decreases since 
this ratio then represents a smaller field as depicted in the 
lower right panel of Fig. [2] However, for the same reason, 
greater magnetic flux ratios become accessible for larger 
rings. 

It is often more instructive to look at the results shown 
in the right panel of Fig. |2] where the familiar diamond- 
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FIG. 3: Plots of the conduction versus magnetic flux for the 
analytical model (upper panel) and for an atomistic calcula- 
tion based on a tight-binding model (lower panel). The ring 
radius in both cases is i? = 5nm. Other parameters in the 
analytical model are s — m* = 4me and Xf ~ 12A. 
Tight- binding parameters are defined in the text. 



shaped conduction pattern within the energy/phase ra- 
tio plane behaves more simply and predictably, in order 
to gain a better intuitive understanding of the system. 
Here, the two ends of the white conduction lines rep- 
resent up/down conduction where the chemical poten- 
tial and magnetic field are at the middle of such lines. 
Varying e still affects only the peak widths, modifying R 
changes the lengths of the conduction lines, changing m* 
scales the diagram in the energy axis without changing 
the length of the conduction lines and altering the energy 
of the conduction electron scales the diagram along with 
the lines in the same axis. Thus, in order to construct an 
efficient spin filter one has to construct this diagram and 
control the position of the white conduction lines by the 
application of a proper magnetic field and gate voltages 
to achieve a desired behavior. 

Good molecular spin filters can be made when two con- 
duction peaks with different spins are sufficiently sepa- 
rated in energy and magnetic field to allow only electrons 
of one spin to traverse the ring. This is best achieved for 
molecular devices with a large effective mass for the con- 
duction electron and small couplings between the leads 
and the ring ("bad contact"). In Fig. |3] we plot the 
magneto-conduction for a molecular device with rela- 



tively high effectiveness. We zoom on the region of the 
first magnetoconductance peak. The upper panel shows 
the results for a analytical model discussed above and the 
lower panel shows the results of an atomistic calculation 
for a molecular AB interferometer as illustrate in the up- 
per panel of Fig. [l] As clearly can be seen, depending 
on the value of the magnetic gate, one can open output 
channel for spin up or down while at the same time close 
this channel for the other spin. 

To calculate the magneto-conductance of a molecular 
AB interferometer we have used a simple tight-binding 
model where we assume a single electron in a spherical 
s-level for each site on the molecular ring. We add the 
proper magnetic terms to the tight-binding Hamiltonian 
of the system: 

= Htb + a*bL ■ B + ^niefilB^Ri + gfinO- ■ B, (12) 

where L is the angular momentum operator, B is the 
magnetic field vector, and Rj^ the projection of R onto 
the plane perpendicular to B. A gauge invariant basis is 
used to evaluate the tight-binding Hamiltonian matrix: 
|l'S)a^ = |ls)ae~Tf -^a-"", where \ls)a is a Is type orbital 
centered on site a, and = — ^(Rq, x Bq) is the vector 
potential evaluated at the position Rq, of site a. We 
take the diagonal matrix elements of Htb to be equal 
to zero (energy scale) and the off-diagonal elements are 
proportional to the overlap between the gauge invariant 
basis on the different electron sites, as described in more 
detail in Refs. [MlandlSTl 

The conductance is calculated using the Landauer for- 
malism^ which relates it to the scattering transmittance 
probability through the system. The transmittance is 
given by T{e) = iTr[GHe)T i{e)G{e)To{e)]. Here, G(e) = 
[eS — Hs + i{Tj + To)] ^ is the retarded Green function, 
S is the overlap matrix, and F/^o are the imaginary parts 
of the self-energy (E) of the input/output channels. For 
the results presen ted in Fig.[3]we use both imaginary ab- 
sorbing potential^2illII and an itera tive semi- infinite bulk 
Green functions calculation schemJ^HinESIlIl to calculate 
the self-energies of the leads. 

Comparing the results of the simple analytical model 
to the results obtained from the tight-binding model in- 
dicates that the same physical picture emerges for the 
tight-binding approach. This is expected based o n previ- 
ous studies where the Zeeman effect was neglected. S^I^IIIZI 
The agreement between the two approaches indicates 
that the diffraction pattern is insensitive to the pertur- 
bations caused by an ionic potential and the results will 
not be invalidated by a more thorough (single-particle) 
treatment. The only free parameter used in the analyti- 
cal theory is the scattering amplitude -^/e, which was ad- 
justed to match the width of the conduction peaks. The 
effective mass entering the analytical model can be cal- 
culated directly from the tight-binding parameters. For 
a single s-level within tight-binding for a ID crystal with 
inversion symmetry and a site distance of a one can show 
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that the dispersion relation can be approximated by^ 
ekO^Eo + a^j{a)k^. (13) 



This can be compared directly with a nearly free electron 



of mass m* and energy = i?o + 5^^^' allowing us to 



identify: 



27 (a) a2 



(14) 



Here 7(a) is the resonance integral, i.e., the off-diagonal 
nearest neighbor matrix element of the tight-binding 
Hamiltonian. 



C. Spin Splitter 

The physics of the three-terminal device shown in the 
lower panel of Fig. [T]was recently discussed for the case 
where the Zeeman splitting was ignored.-- Hod et al. 
showed how cyclic molecular rings can be used as parallel 
magnetoresistance logic gates (in contrast to the switch- 
ing devices based on two-terminal rings ).ESlThe basic idea 
was to couple the cyclic molecular system to three leads 
creating a three terminal device and to apply an external 
magnetic field. By carefully selecting a narrow resonance 
through which conductance occurred they showed that 
such a setup can be used to simultaneously switch one 
channel "on" and at the same time switch the other chan- 
nel "off" . This was achieved by carefully adjusting the 
phase of the conducting electron with diminishing ampli- 
tude on one exit channel and a large amplitude at the 
other channel. A proper combination of a gate potential 
and realistically low (compared to the full AB period) 
magnetic fields were used to obtain parallel logic oper- 
ations such as AND and AND-I-NOT. This was demon- 
strated for a molecular system composed of conjugated 
Benzene rings and further discussed in terms of a single 
channel continuum model. 

The goal of the present study it to develop a spin split- 
ter device by extending the approach presented in Ref . |^ 
to include the Zeeman effect. Specifically, we will show 
how the phase of the conducting electron can be tuned 
such that one exit channel is turned "on" for up spins 
( "off" for the down spins) and at the same time the other 
exit channel is turned "on" for down spins ("off" for up 
spins) . To achieve this we need to increase the parameter 
space necessary to provide means to control the efficiency 
of the device. As will become clear below, the angles be- 
tween the different channels (a, and 7 = 27r — a — /3, 
cf. Fig. [T|) will be used as control parameters to access 
the many ways by which the spin splitter device can be 
implemented. 

We begin with a close examination of the results for the 
three-terminal Zeeman spin splitter described by the sim- 
ple continuum model. The conduction given by Eqs. ( 11 ) 
and ([7])-(|9| for the case of zero bias is plotted in the up- 
per panels of Fig. [4] for the two output channels for the 



case of spinless electrons. Following the analysis of the 
spin filter shown in right panel of Fig. [2] we observe the 
familiar diamond-shaped conduction pattern within the 
energy/phase ratio plane for each output channel. As be- 
fore, the two ends of the white conduction lines represent 
up/down conduction where the chemical potential and 
magnetic field are at the middle of such lines. The ring 
parameters are R = 20nm, m* = 125me, and ^/e = j^. 
We take a = 15.9 deg and /3 = 19.7 deg such that the 
conduction through one channel is related to that of the 
other by a mirror symmetry around the field Bq « 5Tesla: 
gOi{B,e) ~ gOiiBo — B,e). These angles are obtained 
through an optimization procedure to achieve maximal 
effectiveness. 

A spin splitter is obtained when one end of the Zeeman 
split line conducts for one channel only (the up spin for 
the upper left panel of Fig.|4| while at the same time the 
other end of the Zeeman split line conducts for the other 
channel only (the down spin for the upper right panel 
of Fig. [4]). The structure of the conduction allows this 
for points where the Zeeman separation is equal to the 
separation between two ring energy levels that conduct, a 
very stringent requirement. If a satisfying choice of lead 
angles is always available (which appears to be the case), 
then using the simple analytical model of the ring this 
happens when: 



'^3^^^B=,^^[{m + 2f-{m-2f\ (15) 

Or, rewriting B as ^^-^ and /Xg as 1^^, we find that 
the condition where the Zeeman separation is equal to 
the separation between two ring energy levels 



(p 2me 2TOe ttR 

— = mp = , 

Po m*g gm* a 



(16) 



where, as before, a is the separation between the sites on 
the ring. The condition given by Eq. ( 16 1 is a necessary 



but not sufficient condition for a spin splitting device. 
When the condition is met, the angles between the input 
and output channels are optimized to achieve a desired 
efficiency such that a mirror symmetry around the mag- 
netic field given by Eq. ( 16 1 is achieved. 



Similar sets of configurations with odd or half-integral 
flux ratio can be found in much the same way, due to 
the AB splitting. Clearly for rings with nanometer scale 
dimensions, where a flux ratio of one corresponds to 
thousands of Tesla, this is not a viable option unless 
TO* 3> TOg. However, at tens of nanometers many config- 
urations are possible at fields of a few Tesla. The main 
problematic issue that remains is the fact that in order 
to conduct through very low levels on the ring in this 
simple model, the Fermi energy must be lowered very 
significantly. If we assume that the Fermi level is of the 
order of — , one can show that in this simple model the 
magnetic field required to build this device at the Fermi 
energy becomes reasonable only when R is of the order 
of micrometers. However, if one manages to increase the 
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FIG. 4: Transmittance through a three-terminal ring device as a function of the flux ratio (magnetic field) and the kinetic 
energy of the conducting electrons. Left and right panels correspond to output channels 1 and 2, respectively. The ring has 
a radius of J? = 20nm, which corresponds to an AB period of about STesla. Upper panels: m* = 125me, zero bias, and the 
two endpoints of the white line represent electrons of the two spins. Lower panels: m* = lOme, small bias of 6meV, and here, 
spin-independent conduction occurs all along the purple line, while along the white (black) lines, still of length gfis, only up 
(down) electrons are transmitted, as their conduction windows shift in energy according to their spins. Rainbow color scheme 
used (red indicates high transmittance of « 1) . 



spin g-factoi'^ or the effective mass in the device m* (see 



Eq. (14)), a nanometric device would be feasible at mere 
Teslas. This may hint that here as well a realistic de- 
vice is a matter of the choice of materials - one would 
in principle need to custom-design a system in which the 
density of states is high enough at low kinetic energies 
that only levels of such low energies are occupied in the 
ground state, or where conduction sites are far apart and 
weakly connected.'^^'^ 

The problems just mentioned are no longer present if 



one considers biased conduction (we assume the voltage 
falls symmetrically across the junctions to avoid compli- 
cation) as shown in the lower panels of Fig.|4j The appli- 
cation of a finite bias voltage allows realistic configura- 
tions with the same structure at various effective masses. 
The example shown in the lower panels of Fig. [4] corre- 
spond to m* = IOtoe and a bias voltage of 6meV. The 
conduction window can be tuned to contain two quite 
distant levels. If levels with the transmissive properties 
of the ones previously discussed are selected, the Zeeman 
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10.5 11.5 12.5 10.5 11.5 12.5 10.5 11.5 12.5 

B[Tesla] 

FIG. 5: Plots of the spinless current versus magnetic field and gate voltage in the tight-binding approximation for the two exit 
channels of a 20nm 3-terminal ring device (left and middle panels). The spin dependent current is represented by the value of 
the current at the ends of the black lines. Right panel show the corresponding current effectiveness. In the right panel, color 
scale adjusted for maximum detail, such that the red areas represent a maximum effectiveness of approximately 0.6 (rainbow 
color scheme used). The best polarization achieved here is « 90% for one channel and ~ 75% for the other. 



energy need not span the space between them. This can 
be done by the bias voltage and the Zeeman term must 
be no more than the level broadening (or a few k^T if 
this is more) to ensure conduction of only one spin per 
level. However, the conduction of electrons in the entire 
window of bias voltage should vanish in order to make 
efficient spin-splitter devices, since the biased current is 
the integral of the transmission over this energy range. 

In an extension of the graphical method shown, the 
biased conduction window could be represented by a line 
of length eV-B (purple line), where Vb is the bias voltage. 
All "spinless" levels on this line conduct. Centered on 
the line's ends are two parallel lines of length 2g/iB-B 
(black/ white), where only one spin level is shifted into 
the window - therefore, at these ends only one spin type 
takes a part in conduction, while on the rest of the lines 
both do. Now, if this line is placed like the white lines in 
the upper panels of Fig. |4] with the levels at the ends in 
the one spin zones, the desired result is achieved as long 
as the conduction is zero for all other energies covered 
by the bias voltage. This can, however, be done at much 
higher energies and lower fields than those shown in the 
upper panels, since the level spacing spanned by the bias 
voltage can easily be orders of magnitude greater than 

In order to make a stronger argument that the spin- 
splitting configurations are a physical phenomenon rather 
than an idiosyncrasy of the simple one-dimensional an- 
alytical model used here, one might proceed by repro- 
ducing them in a more detailed numerical model. While 
remaining in the independent electron picture, an ob- 
vious and flexible choice is a ring of atoms represented 
by a tight-binding Hamiltonian, for which conductance 
can be calculated with the methods described above. 
The wealth of parameters makes it problematic to fine- 
tune exactly corresponding arrangements between the 



two models, and yet it is quite a simple matter to con- 
struct in one model a ring and lead configuration which 
worked well as a splitter in the other, and to try and see 
how well it works at some choice of parameters which 
should have analogous physical meaning. This should be 
more similar to what might be done in an actual experi- 
ment. 



Since in the analytical model of a 20nm ring we have 
already located, as shown in Fig.|4] a good splitter config- 
uration near lOTesla, the same configuration (with some 
effective mass) would be a convincing place to look for a 
splitter using the tight-binding model. We focus on the 
more realistic low effective mass regime where a bias volt- 
age is needed and on a finite temperature of IK, which 
should still leave the system well within the quantum 
regime. In Fig. [5] we plot the current as a function of the 
experimentally available parameters - the gate voltage 
and the magnetic field strength. The familiar asymmet- 
ric structure from the previous plots is blurred, but still 
readily recognizable in the two left panels, while the right 
panel shows the high effectiveness which can be reached 
when the magnetic field, the gate voltage and the bias 
voltage are all appropriately tuned. The location of the 
effectiveness peak could be easily predicted by the trans- 
mission line method previously discussed. Obviously, it 
would be practically impossible to stumble upon such a 
fortunate combination of conditions by accident, and any 
experiment must search for them under the explicit guid- 
ance of a model such as the analytical model suggested 
above. 
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IV. COULOMB BLOCKADE CALCULATIONS 

A. Charging Energy 

So far, we have ignored the energy it takes to inject 
multiple charges into the small region of the ring. This 
is of some concern since, as will be discussed below, the 
charging energy can be very significant here and espe- 
cially when ring-leads coupling is weak, charging efl^ects 
can play a major role. To stress this point, consider a 
device (ring) that is charged with an electron of spin 
up. The energy it takes to bring another electron with a 
spin down can be significantly different than in the single 
electron picture. Can this additional degree of freedom 
provide means to construct more efficient molecular spin 
filter and spin splitter devices, or will it prevent their 
realization? 



a 
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FIG. 6: The charging energy of a nanometric ring is plotted 
against the ring's geometric parameters. The model used is 
a cylinder of radius b and height h with a cylindrical hole of 
radius a, where b is taken as a + h to provide a square profile. 
Note that the sharp divergent peak at the origin has been 
cropped at i7 = 2eV in order to show greater detail. 

In general, under our set of assumptions one can ex- 
pect that when charging effects are neglected, states on 
the ring will be grouped into pseudo-bands with the same 
total occupation, having widths of around a few level sep- 
arations and separated by the order of the Fermi energy 
(which is usually much greater than the level spacing). 
The addition of the charging term will shift non-neutral 
bands up linearly in the charging energy with an increas- 
ing slope for charged configurations. One should be able 
to compensate for these changes by using a gate or bias 
voltage, as discussed below. 

One of the simplest ways of taking the charging energy 
into consideration in the regime of strong charging energy 
and weak lead bin ding is within a multi-ele ctron master 
equation approach.l^ ^^ l ^^ l ^^ l ^'^ l ^^ l ^m^ IMl The correct- 
ness of the method requires that the effects of broadening 
be negligible, which is quite generally not true, but the 
approximation nevertheless provides some insight into 
the behavior of nanometric devices when charging has 



been accounted for, particularly when a nonzero bias has 
been applied. This formalism is statistical only, and un- 
fortunately this makes discussion of perfect devices that 
send every electron to the right place more difficult than 
it is in the single-electron formalism. Alternatively, one 
can utilize a perturbation treatment in the lead-device 
coupling, recently suggested by Konig and Gefen for an 
AB mesoscopic ring,^^ however, this approach is limited 
to off-resonant transport only. 

Before we proceed to discuss the calculation of the cur- 
rent within the multi-electron master equation approach, 
we brieffy digress to the matter of estimating the charging 
energy and its variation with ring geometry. Calculating 
the energy needed to (uniformly) charge a ring or torus 
is an elementary electrostatics problem, though not one 
to which an analytical solution exists to our knowledge. 
The calculation below provides a qualitative estimate of 
the magnitude of the charging energies and the numbers 
should not be taken as being accurate to within more 
than an order of magnitude. 

The details of the calculation are as follows. As a 
model we used for simplicity a cylinder of radius a and 
height h with a cylindrical hole of radius b through its 
axis (see inset of Fig. [6] for a sketch of the model). The 
potential (j>{r, z) at r and z can be obtained following a 
standard procedure ,1521 and is given by: 




\J {r — r' cos 4>Y + sin^ + (z — z')^ 



The charging energy for a single electron may then be 
found by setting p to ^^^^2_fc2) inside the ring and zero 
otherwise: 

U^ = \j drp0(r, z) (18) 

The charging energy for N electrons is then given by 
Uq{N — NqY where A'o is the neutral number of electrons. 

In Fig. [6] we plot the charging energy of a nanometric 
ring as a function of the cylinder radius a and height h. 
For typical rings considered in this work with a radius of 
a = 5nm, a height oi h= Inm, and h = a + h = 6nm the 
charging energy is of the order of Uq < ^eV. For smaller 
dimensions the charging energy increases considerably 
and can exceed several electron volts for sub-nanometer 
rings. 

B. Calculating the Spin Current 

If one neglects spin-dependent multi-electron effects, 
then it is formally straightforward to construct from the 
set of one-electron Hamiltonian and spin eigenfunctions 
an anti-symmetric basis of multi-electron wavefunctions: 

*nin2...n, =^l...i Vn, (19) 

ni = l 
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Here is the antisymmetrization operator and the 

states are identified by their (spin-dependent) level oc- 
cupations rii (0 or 1 for fermions). Using this anti- 
symmetric multi-electron wavefunction we can uniquely 



and conveniently determine the matrix elements of a 
general many-body operator G according to the Slater- 
Condon rules where only single electron integrals are 
taken into account: 



{iPr\G\ipj) = gij 



M 
j 



gkk 



|G| 

in2...nj ...nk...nM 

\G\ 

^ 711712. ■ .n - ■ . 

{'^ni7i2...nj ...7ik...ni...nM |C|^'„j„2...„' „' = 0, 



(20a) 

(20b) 

(20c) 
(20d) 
(20e) 



where M is the number of single-electron levels taken into 
account, and n'^— \1 — ni\. Multi-electron effects will be 
considered only in the form of charging energy. Since 
these values will be used in a rate-process calculation 
rather than a full quantum formulation, constructing the 
multi-electron states themselves is actually redundant. 



and Eqs. (20bl-(20el along with the single particle data 
will provide all the necessary information. 

In order to perform a master equation based 
estimation of the current, the transfer rates be- 
tween dif ferent m ulti-electron states must first be 

Calculated.™l35l36 37 38 .39 40 41 H ^^^^^ ^^^^ 

only levels near the Fermi energy will take part in conduc- 
tion. This implies that the N levels closest to the Fermi 
energy will be used to construct the 2^ multi-electron 
states themselves using Eq. (20b I, with the additional 
charging term Uq{N — Nq)"^. 

The transfer rate through lead £ between two multi- 
electron states is given by: 



Ri 



r 



e,al3 



(21) 



We label multi-electronic states by the Greek indices a 



(13) 



'AT 



}. Single 



{nl"\^"^---n^^} and /? = {nf,n 
electronic levels are labeled by the indices i and j. We 
also define the total transfer rate summed over all leads: 



Ra^p — ^ Rl,a^f: 



(22) 



In the above equations, F^.q^ is related to the the imagi- 
nary part of the self-energy in the single electron picture. 
To lowest order, F^ = 7^ if the two multi-electronic 
states differ only by the occupation of level z, Tg^afs — Ji.ij 



if they differ only by rii and Uj and Ui 



= 1. Oth- 



erwise F 



e,al3 



= 0. 



rules (cf., Eqs. (20c I and (20d) 



This follows from the Slater-Condon 
^i^ij is the matrix el- 



ement of the imaginary part of the self-energy. Qap in 



Eq. (21) is related to the Fermi-Dirac function, /(e): 

UP^ \l- f{ea-ep- l^e) < N^, (23) 



where Na — J2i "-i"' ^^'^ number of electrons in state 
a. 

With the rates known, the linear master equation sys- 
tem can be written as the condition for steady-state: 

Rc^isPo. - Rp-"-Pp = 0' (24) 

13 f3 

where Pa is the probability that the system is in a multi- 
electron state a. Once the steady-state occupation prob- 
abilities have been s olved f or the current can be expressed 

^ ^33 34 35 36 37 3 8 39 4 01411421 



where 



le = — e Rl,a^i3PaSal3, 



-1 Na<Np, 

Sap = \ -l Na> Nf}, 
Na = Np. 



(25) 



(26) 



Intuitively, this expression states that current flows out 
of lead £ whenever an electron flows from it into the 
device, with the inverse also true. Following a similar 
line of physical reasoning leads to an expression for spin- 
polarized current: up or down current flows out of lead 
£ whenever an up or down electron flows from it into 
the device. Assuming no coupling between levels with 
different spin, the spin-dependent current is given by: 

Ie,m) = ~eY,Ri,a-*pPaST(i)ap, (27) 

ap 
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FIG. 7: Panels (a)-(d): surface plots of the current versus gate and bias voltages. (a),(c) and (b),(d) panels correspond to 
up-spin and down-spin currents, respectively. The lower panels zoom into a region of the voltage plane in which a realization 
of an up spin-filter is possible. Panels (e)-(h): a cut through constant gate voltage plane of the current (e), the effectiveness 
(f), the conduction (g), and the state population (h) as a function of the bias voltage. In all panels the ring is composed of 
40 single electron sites with metallic leads at opposite ends, and a magnetic field of lOTesla is applied. The charging energy is 
^eV and T = |K. 



and 



'Ta)a/3 



+1(0) S^<Sp^N^<N|i, 

0(+l) Sc,> Sfi^N^< Np, 

0(-l) Sa.<Sp^N^> Nfi, 

-1(0) Sc,> SfihNa> Np. 



(28) 



Here, So, — sf^ where si""* — ±1 for spin up or down, 
respectively. 

Peaks in the differential conduction as a function of 
the bias voltage V-b can be expected whenever there ex- 
ists a difference in energy between two states differing 
in their number of electrons by one, which is occupied in 
one lead but not the other, i.e., when the conduction win- 
dow grows to contain a spectral line. This is why levels 
spaced more than about eVe + k^T from the Fermi level 
should not take part in conduction within this formalism: 
transfer through them involves electrons or holes (refer- 
ring here to level vacancies) not present in the leads. 



C. Spin Filter 



A filtering device remains straightforward in this for- 



malism, and the previous discussion in subsection |IIIB| 
pertains to it as well. Several differences are nevertheless 
evident. In the single electron formalism a device will 
conduct at zero bias when the chemical potential coin- 
cides with an energy level on the device. In the master 
equation formalism this still happens when the charging 
energy is zero. Under such conditions, the difference in 
energy between the neutral state and the first charged 
state always equals the Fermi energy, and this difference 



is obviously contained in the zero-bias conduction win- 
dow. 

In the presence of charging, the picture is somewhat 
more involved. The same states discussed above are now 
shifted differently by charging since they have a different 
number of electrons. For the sake of clarity, we will first 
consider the simplest case where only 2 levels and thus 
4 states are included. We assume that the system can 
have A'o — 1 electrons, A'o electrons, or TVg + 1 electrons 
with energies t/o, eig^gB, and 2e + Uo, respectively (e is 
the single electron level energy). Furthermore, Nq is the 
neutral occupation and we assume it is such that one of 
the aforementioned levels is occupied when the systems 
is neutral. The value of these energies was calculated 
by taking the single electron level e and adding charging 
terms {Uq for both Nq — 1 and A^o + 1) ^i^d magnetic field 
splittings (for the A'o state). 

When the chemical potential in one of the leads equals 
the energy difference between multi-electron states a 
peak will appear in the differential conductance. For the 
case described here when the bias voltage falls symmet- 
rically on both junctions, conduction peaks will occur 
when the energy difference between two multi-electron 
states of different occupations is equal to the value of the 
lead Fermi energy combined with half the applied bias 
voltage. Taking the Fermi energy as that of the lower 
of the two single electron conduction levels, the conduc- 
tion peaks are expected to appear at eV^ = 2Uo and 
cVb — 2{Uo + 2g^BB) (these peaks are doubly degener- 
ate). 

Typical behavior of a spin filter device is illustrated in 
Fig. m where Uq = j^eV, e = leV, and T = ^K. On 
the left, surface plots of the current as a function of the 
gate and bias voltages are shown (panels (a)-(d)). Upper 
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panels show the current of spin up (left) and spin down 
(right) electrons for a wide range of gate and bias volt- 
ages. Note how current steps (or conduction peaks) occur 
whenever the conduction window cuts across a difference 
between state energies as discussed above. For instance, 
with no gating transitions with ^ ± 1 occur simul- 
taneously at eVg = 2C/o, and with a gate voltage of ±Uo 
the first step occurs at zero bias. On the scale of the plots 
shown in panels (a) and (b), the currents for the two dif- 
ferent spins are almost indistinguishable. The difference 
in current between the two polarization can be seen in 
panels (c) and (d), where we zoom in on a specific area 
of the fine described by eVc = Uq — 6^3/2, where Vq is 
the gate voltage. Here, one spin starts to fiow at slightly 
lower bias voltages, the current rising with the voltage 
to a high peak and then falling back down as the other 
spin begins to fiow as well. A plane cut through this 
surface is shown in panel (e), where it is clearly demon- 
strated how by tuning the voltages a spin filter can still 
be constructed when charging is taken into account. A 
symmetric filter for the opposite spin type is found when 
inverting the gate voltage. As can be seen in the figure, 
the first conduction peak occurs at eVg = 2{Uo — 2g^^B) 
and the second occurs at eV^ = 2C7o- As a result of the 
broken symmetry, these peaks are non-degenerate and 
correspond to transition with A^ N — 1. The other 
two peaks corresponding to A^ ^ A^ + 1 will occur at a 
significantly higher bias voltage. This is in contrast to the 
simple example where only 4 states were considered, as 
describe above, where the current peaks are degenerate 
and appear simultaneously. 



One of the interesting features shown in Fig. [7] is the 
negative differential spin conductance of the spin up elec- 
tron (panel (g)). This is explained by the sudden drop 
in the population of the |01 > state (corresponding to a 
conduction electron of spin up) as the change in chem- 
ical potentials begins to allow the population of state 
1 10 > (corresponding to a conduction electron of spin 
down). This population switching is reminiscent of the 
nonmonotonic change of occupation in two e lectrostat- 
ically coupled single-level quantum dots.'^SED The state 
populations are shown in panel (h) of Fig. |7] The drop 
in the spin up state population is correlated with the 
occurrence of negative differential spin conductance (the 
conduction is normalized to the maximal valu^^. This 
can be explained in the following way: the current for 
each spin is determined by the product of the probabil- 
ity that the system is in state a and the rate of transitions 
between state a and state |00 > (which is the state with 
no conducting electrons), where a is |10 > for spin up 
conducting electron and |01 > for spin down conducting 
electron. We therefore expect that at chemical potentials 
where the relevant Fermi functions and hence the rates 
are nearly constant, the current will be approximately 
proportional to the population. 



D. Spin Splitter 

The discussion of the spin splitter within the multi- 
electron master equation is more involved than the spin 
filter and requires at least 4 levels (16 states). We as- 
sume that at A'o (the neutral occupation) two of the four 
conduction levels are occupied. Having previously es- 
tablished a way in which single electron levels suitable 
for spin splitting can be found on a ring, we focus on 
the effects of charging. First, and for the moment ne- 
glecting charging energies, we assume that we have found 
two levels with energies e ± A that are spin degenerate 
without Zeeman splitting, with the following properties: 
they should be adjacent (other than perhaps for non- 
transmissive levels, which will not qualitatively affect our 
results), and should both be coupled symmetrically to 
the input lead, while each is coupled to a single different 
output lead. Other levels should be at least 2A away 
in energy. All couplings are taken to be the same and 
diagonal couplings are ignored. Accounting for Zeeman 
splitting will simply shift each of these levels in energy 
by ±g^^B without modifying the coupling. As before 
for the single electron picture, the coupling between lev- 
els with different spins is neglected. An example of such 
a case is discussed above in Sec. |IIIC| where A is given in 
terms of the ring level spacing (c.f., Eq. 15 and Fig. |4]). 



Left Channel 



Right Channel 




FIG. 8: A spin-splitter made from two levels leV apart where 
A — |eV, each of which is further split by the Zeeman energy 
at lOT. The upper diagrams show the conductance through 
the two output channels when the charging energy is ignored 
and the lower diagrams shows the conductance of the two 
output channels at a charging energy of Uo = 0.6eV. We 
have adjusted the gate potential so that in both cases the 
current rises at Vb ~ IVolt and at this value only one spin 
conducts through each lead. 

Setting the chemical potential of the input lead to 
III = e + and both output chemical potentials to 
jiQ = e — it should now be clear under both for- 
malisms that without charging energy a perfect spin split- 
ter exists at eVe = 2A — 2g^^B, as shown in the upper 



14 



panels of Fig. [Sj In the single electron picture this hap- 
pens as two levels with opposite spins and different lead 
bindings enter the conduction window. In the multi- 
electron master equation this happens when the differ- 
ence in energy between N — Nq states and = iVo ± 1 
states enters the conduction window. These two condi- 
tions are equivalent for Ug = 0. 

Including charging effects will cause the = iVo ± 1 
states to move up in energy by Uq. This has the ef- 
fect of repositioning the conductance peaks at eVs — 
2A ± 2gfiBB ± Uq. If Uq is of the order of the Zee- 
man term then the charging energy will interfere with 
the structure that allowed our spin splitter for Uq = 0, 
since the energy difference between states that differ by 
one electron of both spin up and spin down can become 
similar. On the other hand, if the charging energy is in- 
creased further, the transitions N = Nq N + 1 and 
N = Nq ^ N —1 become distinguishable in energy. One 
way to take advantage of this is to apply a gate potential 
in order to shift the state energies by —Uq so that only 
the = iVo ^ -I- 1 transitions are in the conduction 
window (new transitions may appear at similar energies 
if other levels are too close, but since level spacing for a 
ring generally increases with energy we can neglect such 
contributions). Once again at eVs = 2A — 2giiBB perfect 
splitting occurs, though at a smaller total current, since 
fewer state transitions are involved. This is apparent in 
the lower panels of Fig. [8] 



V. CONCLUSIONS AND DIRECTIONS FOR 
FURTHER STUDY 

We have investigated several ways in which rings or 
ring-like structures with a radius of the order of nanome- 
ters, coupled to metallic leads, might be used to construct 
simple spin-sensitive devices. We have focused only 
on the Zeeman splitting to differentiate between spins, 
since it is always present while other effects commonly 
utilized in microscopic structures, such as the Rashba 
and Dresselhaus effects, are generally absent in molec- 
ular/nanometric structures. We believe the niche this 
work occupies in the search for nanospintronic devices is 
yet unexplored: very small coherent structures and weak 
spin-dependent effects have not received much attention, 
despite their formal simplicity and significance. What we 
have shown here is that although building devices under 
the burden of such limitations is difficult, it is possible. 
Considering the scientific and technological benefits of 
such devices we believe it will also be worthwhile. 

The basic calculations we have performed are enough 
to point one in the right direction as to the desired prop- 
erties of molecule-sized Aharonov-Bohm spin devices and 
the conditions which their desired operation might be ob- 
served. A similar methodology can be applied to more 
complex devices or sets of devices. The actual devices 
specifically discussed here were two of the most basic 
- a spin-filter and a spin-splitter or Stern-Gerlach de- 



vice. However, the conclusions drawn and principles laid 
out may easily be extended to many interesting systems, 
from quantum gates (since in theory we can use interfer- 
ence devices to perform general unitary transformations 
between input and output gates) to molecular memory 
(since an electron trapped on the device will modify its 
electrical properties and thus may be detectable at a later 
time) . 

Several complementary methods were utilized during 
the course of this investigation: first, a simple single- 
electron, analytical model in which the parameter space 
can easily be explored, and thus basic intuition about the 
system can be gained. Second, a tight-binding nonequi- 
librium Green's function treatment, which incorporates 
a more realistic physical structure that can be compared 
directly with experimental data, but still assume inde- 
pendent electrons, was applied to similar systems stud- 
ied within the simple analytical model. Finally, a multi- 
electron master equation approach that can be used to 
examine many-particle effects was grafted onto the tight- 
binding results. Here, we choose to focus on electric 
charging, since its effects are energetically dominant. 

Our analysis of the filter, meant to be a test-case, was 
based on the study of the spinless case where it was 
possible to create very narrow conduction peaks near 
zero magnetic field by combining weak device/leads cou- 
plings and a gate v oltage to shift the conduction res- 
onances to B = Q |24|25 | 26 | 27 | gj^^^g ^^le Zccmau spht- 
ting depends linearly on the magnetic field, a spin fil- 
ter always requires a finite magnetic field, at least high 
enough to separate the spin conduction peaks in en- 
ergy by more than k^T. Control over the position of 
the spin-dependent conduction peaks and their widths 
can be achieved by carefully adjusting the kinetic phase 

q^l^ - ^^2m- (e±|i^0,„) and the coupling be- 
tween the device and the leads, respectively. Unlike the 
spinless case where 0^. — ^\/2m*e, the kinetic phase is 
now a function of the magnetic flux itself, and thus, the 
conduction is not a simple periodic function of the pa- 
rameters. High efficiency spin filters are constructed at 
the highest magnetic field possible, where spin-dependent 
effects are strongest, and ffexibility is gained by selecting 
materials or structures with high effective mass (or large 
ring size, which is not desirable). Charging effects do 
not drastically modify this picture since spin selectivity 
depends only the Zeeman term and the application of a 
gate voltage can compensate for the charging energy it- 
self. However, charging leads to a breaking of symmetries 
and as a result to negative differential spin-conduction. 

The physics of the spin splitter device is similar to that 
described above with the added complication that differ- 
entiated control of the spin-dependent wavefunction at 
the different leads is required. Within the single elec- 
tron picture, a spin splitter device may operate when 
two spin-degenerate levels exist such that one transmits 
through one lead only and the other through the other 
lead only. Furthermore, the level separation should equal 
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exactly the Zeeman splitting. Alternatively, the level sep- 
aration can be compensated for by the application of a 
bias voltage, as long as there are no transmissive levels 
between the two spin-degenerate levels. Here, parameter 
space includes the ring size, effective mass, device-leads 
couplings, and the angles between the output leads. In- 
creasing the effective mass, as before, reduces the level 
spacing on the ring, thereby enabling the construction of 
spin splitters at lower magnetic fields and with smaller 
dimensions. Charging breaks the symmetry and can re- 
duce the overall current through the device. However, 
charging does not abolish the general picture and the 
application of a gate potential can be used to overcome 
most of its effects. 

For both prototype devices we have shown how the 
limits of lead-coupling, system geometry and tempera- 
ture at which one might expect to see the desired effects 
can easily be estimated, as well as several ways in which 
one can implement specific behaviors by systematically 
finding parameters at which they occur, either exactly 
or approximately. We have also shown that even within 
the parameters space of a device formed by a single ring 
with two or three leads, nontrivial behaviors with useful 
properties take place. In at least two instances, we have 
made a case that our findings are physical rather than 
a mathematical peculiarity of some model by reproduc- 
ing them under different assumptions and formalisms. 
We have found an interesting example of negative dif- 
ferential conductance for a spin polarized current with 
a simple explanation. Finally, we have described what 
is, to our knowledge, the smallest Stern-Gerlach appara- 
tus ever reasonably conceived of, and one which is fully 
switchable in directionality at constant magnetic field by 
the application of an external electric field. 

Despite their simplicity, the calculations reported here 
provide several novel predictions that seem to be model 
independent. Nevertheless, more refined models are an 
important set of directions to continue along. In par- 
ticular, drawing from the study of crystals and meso- 
scopic systems, we see ways of enhancing the effect of 
Zeeman splitting which require a more elaborate elec- 
tronic structure description. One such is through the 
effects of spin-orbit coupling and the local spin density 



exchange-correlation energy, both of which have been 
used to explain the giant spin (7-factor enhancement that 
has been observed in mesoscopic structures under certain 
conditions.!^ It is quite possible that the same effect can 
be recreated with discrete levels, although this probably 
requires that only a small number of levels be occupied 
(corresponding once again to very low electron densities). 

While we have studied only two- and three-terminal 
devices, with only one injective terminal, it is reasonable 
that four terminal devices will also be of interest as coher- 
ent quantum gates. One of the reasons we have found it 
worthwhile to draw attention to the importance of many- 
particle effects is our hope that in time-dependent cal- 
culations they may be used to create sequential logical 
behavior without sacrificing coherence - for instance, an 
electron may only be able to enter the device through one 
lead, but when it enters it opens up another lead, and the 
infiow of another electron causes both to be discharged 
through a third lead - a sort of sequential AND gate, two 
of which could form a true AND gate. Of course, to con- 
sider useful computation it is necessary to model entire 
networks of such gates, where input electrons enter on 
one set of leads, propagate throughout the network and 
leave it on another set of leads. Such a network presents 
new delocalized challenges unless devices are somehow 
coupled in such a way that no inter-device interference 
takes place. However, even if the devices within the net- 
work are all "imperfect" devices that only approximate 
logic gates, as long as every device has independently ad- 
justable parameters (like a gate voltage) such a network 
forms a fascinating basis for a model of a quantum neu- 
ral networks. These and related issues are still open for 
future study. 
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